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HE: [ B 49] ES EE Scymnus Kugelann 昆虫 主要 捕食 昨 下 .上 曙 虫 等 害虫 ,是 一 类 经 济 上 重要 的 
天 敌 昆 虫 。 目 前 针对 小 毛 靳 虫 属 的 系统 发 育 研究 尚 属 空白 , 亚 属 之 间 的 系统 演化 关系 尚 不 明确 ,为 
了 建立 合理 的 分 类 系统 , 孟 需 对 小 毛 黑 虫 属 的 亲缘 关系 进行 研究 和 探讨 。【 方 法 】 以 华南 农业 大 学 
4 35,0 ТАЛАҚ 55 亚 属 共 44 种 为 研究 对 象 , 采 用 PCR 技术 对 12S, 165 和 28S rRNA 基因 的 部 分 
序列 进行 扩 增 ;运用 MEGA 7.0 ЖТА RA 125, 165 和 28S rRNA 基因 的 碱 基 组 成 ,基于 
K2P 模型 计算 了 小 毛 款 下 属 44 种 的 种 间 遗 传 距离 ;采用 最 大 似 然 法 (maximum-likelihood，ML ) 和 
贝 叶 斯 推断 法 (Bayesian-inference ВІ) 构建 该 属 的 系统 发 育 树 。【 结 果 】 扩 增 获得 小 毛 默 虫 属 44 种 
的 12S rRNA 基因 序列 平均 长 度 为 356 bp, 16S rRNA 基因 序列 平均 长 度 为 351 bp, 285 rRNA 基因 
序列 平均 长 度 为 315 bp; 序 列 分 析 表 明 ,12S rRNA 基因 的 A, T, G 和 (fC 平 均 含 量 分 别 为 38.8% ， 
43.5%, 11.9% fe 5. 806 , 16S rRNA 基因 的 A,T,G 和 fC 平 均 含 量 分 别 为 37.6% 40.3%, 
14.4% 和 7.7% , 288 rRNA 基因 的 A, T, G 和 fC 平 均 含 量 分 别 为 26.7% , 18.3%, 31. 496 和 
23.596 ;基于 联合 序列 分 析 的 种 间 遗 传 距离 为 0.004 ~ 0.276 ,平均 遗传 距离 为 0.115。 系 统 发 育 分 
析 结 果 表 明 DEI E 2E AR ARR md E SK x TE Scymnus ( Scymnus) Kugelann, ЛАЖ 2 5 
Scymnus ( Neopullus ) Sasaji, АЖ © Æ Scymnus ( Pullus) Mulsant je JA «] S& Ж E 5 Scymnus 
( Parapullus) Yang 均 为 并 系 起 源 。 【结论 ] 基 于 12S, 165 4e 285 rRNA 基因 序列 的 小 毛 昧 虫 属 系统 
发 育 分 析 显 示 传 统 的 形态 学 分 类 体系 与 基于 分 子 数据 分 析 的 结果 部 分 不 一 致 ,这 表明 应 该 对 该 属 
内 各 亚 属 的 鉴别 特征 进行 全 面 检视 ,筛选 并 确立 各 亚 属 的 形态 指标 ,同时 也 表明 该 属 内 的 亚 属 分 类 
单元 需 重 新 厘定 。 
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(Coleoptera ) Z В MV. Н ( Polyphaga ) Д FH. f 


Abstract: [Aim] The genus Scymnus Kugelann is one of the most economically significant genera of 
natural enemy insects due to its predatory activity against pests including scale insects and aphids. 
However, the classification of this genus is still ambiguous, and the relationship of the subgenera within 
this genus is still in dispute. In order to establish a reasonable classification system of the genus 
Scymnus, the phylogenetic analysis is needed to further refine the taxonomic status of the subgenera 
within this genus. [Methods] The partial sequences of 125, 16S and 285 rRNA genes from 44 species 
belonging to five subgenera within the genus Scymnus preserved in the Collection of South China 
Agricultural University were amplified using PCR. The inter-species genetic distance among the Scymnus 
species was calculated based on Kimura-2-parameter model using MEGA 7.0. The base composition of 
128 rRNA, 165 rRNA and 285 rRNA genes was analyzed by using MEGA 7.0. The phylogeny trees were 
generated by maximum-likelihood ( ML) and Bayesian-inference ( ВІ) methods. [Results] The mean 
length of the 125, 165 and 285 rRNA gene sequences amplified from the 44 species of Scymnus are 356, 
351 and 315 bp, respectively. The average contents of A, T, G and C of 128 rRNA gene sequence are 
38.896, 43. 596 , 11. 996 and 5. 896 , respectively, those of 168 rRNA gene sequence are 37. 6%, 
40.3%, 14. 496 , and 7. 796 , respectively, and those of 28S rRNA gene sequence are 26.7% , 18. 396 , 
31.4% and 23. 5% , respectively. The inter-species genetic distance is 0. 004 — 0. 276 based on the 
concatenated sequence data, and the average genetic distance is O. 115. The phylogenetic analysis 
revealed that the genus Scymnus is monophyletic, while the subgenera Scymnus (Scymnus) Kugelann, 
Scymnus ( Neopullus ) Sasaji, Scymnus ( Pullus) Mulsant and Scymnus ( Parapullus) Yang are 
paraphyletic. [Conclusion] The phylogenetic analysis of Scymnus based on 12S, 16S and 28S rRNA 
gene sequences illustrated that the molecular classification system of this genus is not consistent well with 
its morphological classification system. In order to establish a reasonable classification system, other valid 
morphological characteristics should be explored to distinguish the relationships among the subgenera. 
Moreover, more research is needed to further refine the taxonomic status of the subgenera. 


Key words: Scymnus; Coccinellidae; 125 rRNA; 165 rRNA; 285 rRNA; phylogeny 











ЗЕН JE Scymnus Kugelann ЖЕН ШІ 17 5 5. (Neopullus) Sasaji , 0L] NC ХР) S. 








( Parapullus) Yang, S. ( Mimopullus) Fürsch , 直 小 毛 


( Coccinelloidea) B4 E F} ( Coccinellidae ) `6 9 E Wl. BK "p WV JE 5. ( Orthoscymnus ) Canepari 和 S. 
Ж. ( Scymninae) /| E BUE JA ( Seymnini ) ,是 世界 性 分 ( Canalipullus) Lafer( Chen et al., 2013) 。 我 


布 的 类 和 群 。 目 前 全 世界 球 虫 科 种 类 已 知 超 过 6 000 ” 共 记 录 其 中 5 MER , 2r 317] NES rh Ч Je. SUE 
种 (Slipifiski, 2007) ,其 中 小 毛 际 虫 属 已 记录 800 种 亚 属 、 毛 束 虫 亚 属 、 拟 小 球 虫 亚 属 ЛХ BUE IJE 



































以 上 (Chen et al., 2013) 。 我 国 幅员 辽阔 ,地 跨 东 洋 (Chen et al., 2013, 2016) , 

区 和 古 北 区 两 大 动物 地 理 区 系 , 球 虫 资源 丰富 ,已 记 我 国 小 毛 球 虫 的 分 类 研究 起 步 较 晚 , 但 发 
录 1 016 fp, Hirt ES Bm 269 种 , 占 我 国标 虫 科 ，” 速 。 有 学 者 分 别 记 录 小 毛 标 虫族 3 属 4 种 (小 
总 数 的 26.5% ,是 我 国际 虫 科 中 种 类 最 多 的 一 个 属 BRA 2 种 ) ( 刘 崇 乐 , 1963) 8 属 27 种 (其 中 小 毛 


( 庞 虹 等 ，2004; Chen et al., 2012, 2013, 2014, HUR JE 17 Ж) (Pang and Gordon, 1986) ,此 后 
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2015а, 2015b, 2015c, 2016) 。 该 属 已 知 种 类 均 为 。 续 增 加 143 种 成 员 ( 庞 虹 等 , 2004; 虞 国 跃 , 2004) , 
捕食 性 ,主要 捕食 蚜虫 、 内 虫 等 农林 害虫 ,是 当前 我 。 这 些 研 究 工作 为 我 国 该 属 的 系统 分 类 奠定 了 
国 发 展 “ 环 境 友 好 、 生 态 安全 的 害虫 防 控 策略 和 技 ”基础 。 近 年来, 我国 小 毛 际 虫 属 的 分 类 研究 已 取得 


术 "中 害虫 自然 控制 因子 的 重要 力量 。 
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在 世界 范围 内 ,小毛 球 虫 属 已 记述 8 个 亚 属 ,分 (0 Chen et al., 2012, 2013, 2014, 2015a, 20 





ЈА 90 IE JE Seymnus (Seymnus) Kugelann ,小 2015c, 2016), ЖАП, Н FNE BLUE% АХ 
虫 亚 属 S. (Pullus) Mulsant, S. (Поп) Casey, 形态 极其 相似 ,种 间 形 态 差 异 极 小 ,难以 进行 比 





要 的 


较 大 的 进步 ,一 大 批 新 种 和 新 记录 种 陆续 被 描述 、 发 
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究 ,为 种 类 的 鉴定 带 来 较 大 困难 ,从 而 导致 该 属 存在 
许多 近似 种 和 存疑 种 。 尤 其 是 早期 记述 的 种 类 , 描 
述 简单 ,缺乏 雄性 外 生殖 需 研 究 ,部 分 种 类 被 作为 新 
种 重复 发 表 , 一 些 有 疑问 的 同 物 异 名 和 异物 同名 等 
问题 一 直 悬 而 未 决 ( 广 国 跃 和 庞 雄 飞 ，1992; 庞 虹 
等 , 2004) 。 此 外 ,小 毛 球 虫 属 内 亚 属 的 主要 鉴别 特 
征 为 第 1 腹 板 后 基线 完整 或 不 完整 ,触角 101756 
11 节 , 因 鉴别 特征 单一 或 形态 指标 模糊 ,从 而 导致 
亚 属 间 的 鉴别 特征 严重 交叉 重合 ,造成 部 分 种 类 的 
分 类 地 位 混乱 , 亚 属 分 类 阶 元 存在 不 合理 性 及 系统 
发 育 关系 不 明确 , 亚 需 尽快 澄清 。 
线粒体 基因 为 严格 的 母系 遗传 , 因 其 结构 简单 、 
进化 速率 较 低 等 特点 ,普遍 被 用 于 解决 物种 的 起 源 
与 进化 关系 等 方面 的 研究 ( 朱 保 健 等 , 2010). 165 
rRNA 是 线粒体 基因 组 中 既 含 有 高 度 保守 的 序列 区 
域 ,又 含有 中 度 保守 和 高 度 变化 的 序列 区 域 的 一 段 
序列 ,因此 被 认为 是 科 级 以 下 水 平 物种 间 亲 缘 关 系 
和 系统 演化 研究 的 有 效 标记 (Foighil et al., 1995; 
Schneider апа Foighil, 1999) 。 在 叶 蜂 科 ( 王 德 明 ， 
2010), ААЖ Cf ER TÉ, 2009), RER OBW, 
2010) ERPF ( Robertson et al., 2013) 等 昆虫 的 系 
统 发 育 研 究 中 得 到 了 广泛 应 用 。12S rRNA 是 线 粒 
Ж ОМА 上 另 一 个 rRNA 基因 , 因 其 易 扩 增 且 进化 速 
度 快 等 特点 ,也 被 广泛 应 用 于 研究 生物 种 属 间 系 统 
发 育 关系 ( 罗 育 发 和 颜 享 梅 , 2006 ) 285 rRNA Ж 
源 于 核 基 因 ,避免 了 线粒体 DNA. 由 于 单 性 遗传 对 分 
文 分 析 造 成 的 影响 ,也 可 用 于 进行 属 内 不 同 种 之 间 
的 系统 发 育 分 析 (Ohnishi et al., 2003; Cepeda et 
al., 2012) 。 以 上 3 种 分 子 标记 已 广泛 应 用 于 肚 虫 
科 的 系统 发 育 关 系 分 析 中 (Giorgi et al., 2009; 
Magro et al., 2010; Seago et al., 2011; Szawaryn et 
al., 2015) , 
BEFRA HI Ay TAE SIL Жыға АРЫ 
ВЕЛАТ НИЕ T TERR, ЕНЕҢ 
















































































Hi ( Epilachnini ) BUE J% ( Coccinellini ) #125 2 91 
虫族 ( Chilocorini ) 的 系统 发 育 关 系 有 了 深刻 的 认识 
(Szawaryn et al., 2015; Escalona et al., 2017; 李 文 
A, 2017). [He , УЫ SIUE B BS RRA FE EAE 
目前 尚 属 空白 , 亚 属 之 间 的 系统 演化 关系 尚 不 明确 ， 
为 了 建立 合理 的 分 类 系统 , 咀 需 对 小 毛 球 虫 属 亲缘 
关系 进行 研究 和 探讨 。 本 研究 以 华南 农业 大 学 的 馆 
藏 标本 为 材料 , 扩 增 了 小 毛 球 虫 属 5 个 亚 属 共 44 种 
Й] 125, 165 和 28S rRNA 基因 序列 部 分 片段 ,并 且 
从 GenBank 收集 了 其 他 球 虫 和 外 和 群 的 相关 序列 ,分 
析 了 小 毛 球 虫 属 各 亚 属 之 间 的 系统 发 育 关 系 , 其 结 
果 对 厘 订 我 国 小 毛 际 虫 属 内 亚 属 的 划分 等 分 类 问题 
提供 了 依据 ,同时 也 为 建立 更 加 完善 的 球 虫 科 分 类 
系统 提供 了 重要 资料 。 




































































1 材料 与 方法 


1.1 iX 

本 实验 所 用 标本 均 来 自 于 华南 农业 大 学 生物 防 
治 教育 部 工程 研究 中 心 标 本 室 馆 藏 标本 ,标本 保存 
于 无 水 乙醇 中 ,具体 的 采集 信息 及 GenBank 登录 号 
见 表 1。 共 选取 了 小 毛 飘 虫 属 5 亚 属 44 种。 同时， 
选取 了 伪 球 虫 科 ( Endomychidae ) 内 3 个 种 为 外 群 ， 
已 有 的 系统 学 研究 证 明 该 科 是 球 虫 科 的 姐妹 群 
( Robertson et al., 2008; Marvaldi et al., 2009; Zhang 
et al., 2018) 。 此 外 ,还 选取 了 食 螨 球 虫 属 Stethorus 
(1) RISE Е 906 3865 Nephus(3 种 ) 为 外 群 。 由 于 
在 Giorgi & (2009 ) 和 Robertson 等 (2015 ) 的 研究 结 
果 表 明 小 毛 束 虫 属 与 食 螨 际 虫 属 互 为 姐妹 群 , 而 
Magro 等 (2010 ) 的 研究 结果 却 显 示 小 毛 标 虫 属 与 这 
叶 毛 杜 虫 属 互 为 姐妹 群 ,因此 本 研究 同时 加 入 食 螨 
球 虫 属 和 弯 叶 毛 球 虫 属 ,来 验证 其 与 小 毛 球 虫 属 之 
间 的 关系 。 共 51 个 样本 ( 表 1) 用 于 最 终 的 系统 发 


育 分 析 。 




































































表 1 样本 信息 及 GenBank 登录 号 
Table 1 Sample information and GenBank ID 






























































编号 种 名 采集 地 点 经 纬度 GenBank 登录 号 GenBank ID 

No. Species Collecting locality Latitude and longitude 125 rRNA 165 rRNA 285 rRNA 
1 191 Scymnus (Neopullus) fuscatus ”福建 明 溪 Mingxi, Fujian 26°34'№, 116°56'E МК075886 МН742332 МК075842 
2 Ж S. (Л. ) hoffmanni 安徽 石 台 Shitai, Anhui 30960'N, 117°30'Е МК075887 МН742333 МК075843 
3 СЕД 5. (М. ) lijiangensis 云南 昆明 Kunming, Yunnan 26°04'№, 10225076 МК075888 МН742334 МК075844 
4 ЛЕВУ 5. (М. ) minisculus 湖南 炎 陵 Yanling, Hunan 26?30'N, 114°02Е MK075889 МН742335 МК075845 
5 Ж S. (М. ) paralleus 天 津 敬 县 Jixian, Tianjin 40°10'№, 117932  MK075890 МК075930 МК075846 

ЧИ EC XE SUR Taikan Р 

6 西藏 亚 东 Yadong，Xizang 27?25'N, 88°56'E МК075891 МК075931 МК075847 





S. (Orthoscymnus) rhododendri 
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续 表 1 Table 1 continued 




































































































































































































































































编号 种 名 采集 地 点 经 纬度 GenBank 登录 号 GenBank ID 
No. Species Collecting locality Latitude and longitude 125 rRNA 165 rRNA 285 rRNA 
ШІН ШАЙХ: у 
7 S. (Parapullus) alishanensis 台湾 阿里 山 Alishan, Taiwan 23°31'N, 120°48'E MK075892 МК075932 МК075848 
8 Б/Н S. (Ра. ) nanlingicus m 25?48'N, 110°24'Е  MK075893 - МК075849 
9 ЗДАН S. (Ра. ) secula 台湾 阿里 山 Alishan, Taiwan — 23?31'N, 1202487Е  MK075894  MK075933  MK075850 
10 — Bim MALE S. (Pullus) accamptus 云南 景 东 Jingdong, Yunnan | 24?33'N, 101°55 Е | MK075895  MH742336 МК075851 
11 Л 5. (P. ) ambonoidea 云南 景 洪 Jinghong, Yunnan — 21?55'N, 101°09 Е — MK075896 МК075934 МК075852 
12 12/11 S. (Р. ) bourdilloni 8:28 Motuo, Xizang 29°26'N, 95?40'E МК075897 МН742337 МК075853 
13 KENILE 5. (P. ) eryphaconicus 陕西 汉中 Hanzhong, Shaanxi 33°33'№, 107?55'E — MK075898 МН742338  MK075854 
14 ЖЭЛ 5. (P. ) facetus (LASS Мошо, Xizang 29?19'N, 95?20'E MK075899  MK075935 МК075855 
15  *AErRIL/ NBC S. (P. ) fanjingicus 台湾 阿里 山 Alishan, Taiwan — 23?31'N, 1202487Е  MK075900 МН742339 МК075856 
16 Ж/Д: S. (P. ) fujianensis 云南 独 龙 江 27?A1'N, 98°18'°ЕЫ МК075901 МН742340 МК075857 
Dulongjiang，Yunnan 
17 УЖАМ 5. (Р. ) godavariensis 西藏 林芝 Linzhi，Xizang 29?34'N, 95?23'E MK075902 MH742341 МК075858 
18 MRAUK S. (P. ) hingstoni 西藏 林芝 Linzhi, Xizang 29?34'N, 94?29'E MK075903 | MK075936 МК075859 
19 НЖЯ/ЛШЦЫ 5. (P. ) japonicus 云南 动 海 Menghai, Yunnan — 21?35'N, 100°25'Е МК075904  MH742342  MK075860 
20 BRILE S. (Р. ) martensi 西藏 错 那 Cuona, Xizang 27?49'N, 91?45'E MK075905 МН742343  MK075861 
21 ФЛУД S. (P.) nephrospilus 安徽 石 台 Shitai, Anhui 30?60'N, 117°30'E МК075906 МН742344 МК075862 
22 НОЛ ДЈ S. (Р. ) oestocraerus (LASS Motuo, Xizang 29?34'N, 95?23'E MK075907 МН742345  MK075863 
23  5E3L/NSUR S. (Р. ) penicilliformis 山西 垣 曲 Yuanqu, Shanxi 3592077, 112°01'E МК075908 МК075937 МК075864 
24 ШІМ 5. (Р. ) perdere 安徽 黄山 Huangshan, Anhui — 30?40'N, 118°11' Е МК075909 МН742346 МК075865 
25 ДЕЛ S. (Р. ) podoides 安徽 石 台 Shidai，Anhui 30?60'N, 117°30'E МК075910 МН742347 МК075866 
26 ЖЕЛ S. (P. ) posticalis 云南 红河 Honghe, Yunnan 22?53'N, 102°56Е МК075911  MH742348 МК075867 
27 ЖЫҢУМШҢ S. (P. ) scalpratus 台湾 屏 东 Pingdong, Taiwan — 22?12'N, 120°47 Е МК075912 MH742349  MK075868 
28 MNILE S. (P. ) scrobiculatus 云南 动 海 Menghai, Yunnan — 21?35'N, 1002255Е МК075913  MH742350  MK075869 
29 SMALE S. (Р. ) shirozui 云南 动 海 Menghai, Yunnan 22?11'N, 100?20'E МК075914  MH742351  MK075870 
30 xS S. (P.) sodalis 云南 红河 Honghe, Yunnan 22?40'N, 103956 Е МК075915  MH742352 МК075871 
31 BWILE S. (Р. ) takasago 台湾 茂 林 Maolin, Taiwan 22?54'N, 1209407: МК075916 МН742353 МК075872 
32 ФМД S. (P. ) tegminalis 云南 动 海 Menghai, Yunnan 21935"Х, 100?25'E МК075917  MK075938  MK075873 
33  Ng££/j-ESK d 5. (5сутив) acidotus . 去 南西 双 版 纳 22?02'N, 101°01'E МК075918 МН742354 МК075874 
Xishuangbanna, Yunnan 
34 КЕЛ S. (5. ) crinitus 云南 大 理 Dali, Yunnan 25?57'N, 100°23'E МК075919 MH742355 МК075875 
35 ХХА 5. (5. ) dipterygicus 云南 昆明 Kunming, Yunnan — 24?57'N, 102°37'E  MK075920 МН742356  MK075876 
36 ЖЛЛЧЫШЫҢ S. (S. ) dolichonychus 福建 明 溪 Mingxi，Fujian 26?34'N, 116°56'E MKO75921 MH742357 МК075877 
37 INESE S. (S. ) elagatisophyllus — zxW xi Jinghong, Yunnan 21°55'№, 101°09'E МК075922  MK075939 МК075878 
38 КЕЛЛЕ S. (S. ) folchinii 安徽 黄山 Huangshan, Anhui — 30?40'N, 118°11'E МК075923  MK075940 | MK075879 
39 ЖЕЛ S. (S.) inderihensis 山西 宁 武 Ningwu，Shanxi 389507Х, 112°05'E МК075924 МН742358 МК075880 
40 РЖ/ЛМЫШЫҢ S. (S. ) kabakovi B 2. pm 27?Al'N, 98°16'Е MK075925 МК075941 МК075881 
41 VELE S. (5. ) longmenicus 福建 天 宝 岩 Tianbaoyan, Fujian 25?56'N, 117°32'E МК075926 МН742359  MK075882 
42 Х/ЛЧЫШЫҢ 5. (8.) nubilus 台湾 屏 东 Pingdong, Taiwan | 22?12'N, 120°47 Е МК075927 МК075942  MK075883 
43  3EJIVh-ESUR 5. (5. ) scapanulus 山西 垣 曲 Yuanqu, Shanxi 35?20'N, 112°01'E МК075928 МК075943 МК075884 
44 ЖЙ/УЗЫНЫҢ S. (5. ) sufflavus 云南 昆明 Kunming, Yunnan | 26?04'N, 102°50'E — MK075929 МК075944 МК075885 
45 Мерһиѕ includens GU073810 | GU073867 С0073760 
46 Мерһиѕ reunioni GU073811 | GU073868 С0073761 
47 Мерһиѕ bisignatus GU073809  GU073866 С0073759 
48 RARAUHE Stethorus punctillum С(073817 С0073874 С0073767 
49 Trochoideus goudoti GQ302080 С0302160 КР829338 
50 Trochoideus desjardinsi GQ302086 (С0302159 С0302294 
51 Мусеіаеа subterranea GQ302090 С0302161 С0302299 
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1.2 基因 组 DNA 提取 

血液/ 组 织 / 细 胞 基因 组 织 提 取 试 剂 盒 
(TIANamp Genomic DNA Ки) 用 于 总 DNA 提取 ,由 
于 成 虫 体 型 较 小 ,因此 选择 成 虫 的 头 部 或 胸部 肌肉 
组 织 用 于 DNA 提取 ,详细 的 提取 步骤 严格 按照 试剂 
盒 说 明 书 进行 。 提 取 的 DNA 置 于 -20% 保存 备用 。 
1.3 PCR 扩 增 和 测序 

PCR 反应 扩 增 引物 参考 Whiting 等 (1997 ) 和 
Magro 等 (2010 ) 的 设计 ,12S rRNA 基因 扩 增 引物 
125 ai ( 5'-АААСТАССАТТАСАТАСССТАТТАТ-3') 
41125 bi( 5'-AAGAGCGACGGGCGATGTGT-3') ;16S 
rRNA 基因 扩 增 引物 16S F(5'-ATGTCTTTTTGAKWA 
TAATWTAA-3'),168 Б ( 5'-ACGCTGTTATCCCTAA 
GGTAATTT-3^) ;28S rRNA 基因 扩 增 引物 28S sa( 5'- 
GACCCGTCTTGAAACACGGA-3') #1 285 sb(5'-TCG 
САСССААССАССТАСТА-3”) ( Whiting et al., 1997; 
Magro et al., 2010) 。 扩 增 体系 (25 pL) :去 离子 水 
10 uL, PCR Mix 12 р, 上 下 游 引物 (10 umol/L) 4% 
1 uL, DNA 模板 1 kL。 反 应 程序 :94Y 预 变 性 3 
min; 94C 变性 30 s,，50% 退火 30 в, 72% 延伸 1 
min, 35 个 循环 ;72% 最 终 延 伸 5 min, 

取 PCR 最 终 产 物 5 uL 于 绿 如 蓝 核 酸 染 料 
(DNA Green 10000 x ) 染色 的 1% 琼脂 糖 凝 胶 中 电 
泳 ,然后 在 紫外 灯 照 射 下 观察 是 否 有 条 带 , 将 有 条 带 
的 样品 送 上 海 美 吉 生物 医药 科技 有 限 公 司 进行 序列 
测定 。 

1.4 序列 分 析 

序列 的 编辑 、 比 对 工作 在 Genious 7. 1.4( Kearse 
et al., 2012) 软件 中 进行 ,在 拼接 过 程 中 利用 基因 峰 
图 判断 序列 的 质量 ,并 进行 人 工 校 准 , 通 过 该 软件 中 
Trim 功能 切除 两 段 引物 序列 。 使 用 SequenceMatrix 
( Vaidya et al., 2011) 将 3 个 基因 片段 拼接 起 来 , 输 
出 文件 保存 为 FASTA/PHYLIP/NEXUS 格式 进行 后 
续 分 析 。 

利用 MEGA 7.0(Kumar et al., 2016) 软件 对 各 
基因 序列 组 成 进行 统计 ,内 容 包括 变异 位 点 
( variable sites) .保守 位 点 (conserved sites) .简约 信 
































息 位 点 (parsimony information sites)、 自 毅 位 点 
(singleton sites ) 、 碱 基 组 成 、 转 换 与 颠 换 比值 (Ts/ 
Tv) ,然后 基于 联合 序列 使 用 Kimura-2-parameter 
( K2P) 模 型 计算 种 间 遗 传 距离 等 (Kimura, 1980)。 
1.5 系统 发 育 分 析 

首先 基于 单 基 因 分 别 构建 系统 发 育 树 , 碱 基 替 
换 模 型 的 选择 采用 MrModeltest 3.7 ( Nylander, 








2004) 软件 下 的 Akaike Information Criterion ( AIC ) 
评估 选择 最 优 模 型 。 然 后 基于 联合 基因 构建 系统 发 
育 树 ,在 联合 基因 之 前 需要 对 联合 基因 进行 同 质 性 
检验 (incongruence length difference test, ILD ) ( Farris 
et al., 1994) 以 检测 联合 序列 之 间 的 系统 学 信息 是 
否 一 致 ,是 否 适 合 基因 序列 的 联合 分 析 。 检 测 结果 
中 当 P >0.01, 说 明 差异 不 显著 ,可 进行 联合 序列 分 
析 ; 如 果 P<0.01, 则 说 明 序列 不 适合 联合 分 析 。 同 
质 性 检验 使 用 PAUP 4.0 软件 (Swofford, 2001)。 然 
后 在 对 联合 序列 构建 系统 发 育 树 前 ,采用 了 基因 的 
分 区 策略 ,并 与 PartitionFinder 1.1.1(Lanfear et al., 
2012) 中 基于 AIC 准则 启用 全 部 式 搜寻 (All Search) 
寻找 最 合适 的 分 区 策略 和 碱 基 蔡 换 模型 。 

最 大 似 然 树 和 贝 叶 斯 推断 树 的 构建 分 别 采用 
RAxML 8.2.8(Stamatakis, 2014) 和 MrBayes 3. 2. 6 
(Ronquist and Huelsenbeck, 2003) 软件 。 在 利用 
RAxML 8. 2. 8 分 析 的 过 程 中 ,采用 快速 朴 山 法 
( rapid hill-climbing) 搜寻 最 大 似 然 树 ,搜寻 次 数 为 
100 次 。 自 举 值 法 (bootstrap value, ВР) Æ% 1 000 
次 用 于 检验 各 节点 的 可 信和 度 。 在 MrBayes 3.2.6 的 
建树 过 程 中 ,每 次 分 析 都 进行 两 次 独立 的 运算 ,运行 
时 设置 随机 起 始 树 ,每 个 运算 都 包括 4 条 马尔 科 夫 
BEC 条 热 链 1 条 冷 链 ) 同时 和 运行 ,迭代 运算 50 x 10" 
代 , 每 5 000 代 记录 一 次 样本 参数 和 系统 树 。 两 个 独 
立 运行 所 得 到 的 树 中 ,Burn-in 值 设 置 为 25% ,剩余 
的 样本 树 合并 为 一 致 树 。 后 验 概率 ( posterior 
probability, PP) 用 来 检验 贝 叶 斯 推断 一 致 树 各 节点 
的 可 信和 度 。 然 后 在 Tracer 1. 6 ( Rambaut et al., 
2014) 中 评估 马尔 可 夫 链 (Markov chain ) 各 参数 的 
收敛 程度 及 贝 叶 斯 推断 树 的 拓扑 结构 的 稳定 性 。 最 
后 运用 EvolView ( Zhang et а!., 2012; He et al., 
2016 ) 对 系统 发 育 树 进行 信息 标注 。 


































































































2 结果 


2.1 ЛЕД Е 125, 165 和 28S rRNA 基因 序列 
分 析 

共 获 得 131 条 rRNA 基因 序列 ,其 中 12S rRNA 
44 条 、16S rRNA 43 条 和 28S rRNA 44 条 ,只 有 南 岭 
WILE S. (Ра. ) nanlingicus 的 16S rRNA 基因 未 
扩 增 出 来 ,说 明 所 选 引物 具有 很 高 的 扩 增 效率 。 所 
有 序列 利用 NCBI ( https: // www. ncbi. nlm. nih. 
gow) 的 BLAST 进行 相似 性 检索 ,以 确定 所 扩 增 片 
段 为 目的 序列 。 随 后 将 所 有 序列 上 传 至 GenBank, 
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详细 的 登录 号 见 表 1 。 

利用 MEGA 7.0 统计 目的 序列 间 碱 基 转 换 与 颠 
换 的 比值 .差异 百分数 以 及 各 位 点 以 及 A, T, GA 
C 的 平均 含量 。 其 中 12S rRNA 基因 扩 增 长 度 约 为 
356 bp ,保守 位 点 177 个 ,变异 位 点 172 个 ,简约 信 
息 位 点 130 个 , 自 裔 位 点 42 个 ;A, T, G 和 5C 的 平 
均 含 量 分 别 为 38.8% , 43.5%, 11.9% 和 5.8% ,其 
中 Ar+T 人 含量 (82. 3%) 明显 高 于 G+C 含量 
(17.796) , Ez AT Н ЖА/ He (Тв/Ту) 比值 为 
0.517 ,替换 以 颠 换 为 主 。16S rRNA 基因 扩 增 长 度 
2] 351 bp 左右 ,保守 位 点 185 个 ,变异 位 点 171 
个 ,简约 信息 位 点 136 个 , 自 裔 位 点 35 个 ;A, T, G 
和 C 的 平均 含量 分 别 为 37.6% , 40.396 , 14. 496 和 
7.7% ,其 中 A+T 含 量 (77.9% ) 明显 高 于 G+C 含 
量 (22. 196) , 核 苷 酸 的 转换 / 颠 换 (Ts/Tv ) 比值 为 
0. 279 ,替换 以 颠 换 为 主 。28S rRNA 基因 扩 增 长 度 
约 为 315 bp 左右 ,保守 位 点 199 个 ,变异 位 点 115 
个 ,简约 信息 位 点 96 个 , 自 裔 位 点 18 个 ;A, T, G 
和 C 的 平均 含量 分 别 为 26.7% ,18.3% , 31.4% 和 
23.5% ,其 中 A +T & E (45.0%) ЕЕС +С 5 
量 (54. 996) , АМЕ H/E d ( Ts/Tv ) 比值 为 
1. 877 ,替换 以 转换 为 主 。 

基于 12S, 165 和 28S rRNA 联合 基因 序列 使 用 
Kimura-2-parameter (КОР) 模型 计算 了 小 毛 球 虫 属 
各 种 间 的 遗传 距离 ( 表 2) ,其 中 长 毛 小 毛 标 虫 5. 
(S. ) crinitus 55 ХХ / EALE S. (S. ) dipterygicus 
遗传 距离 最 小 仅 为 0. 004 , 南 岭 拟 小 际 虫 S. (Pa. ) 
nanlingicus TED E^ NS S. (Р. ) hingstoni 遗传 距 
离 最 大 为 0.276 ,平均 遗传 距离 为 0.115。 
2.2 基于 12S, 16S 和 28S rRNA 基因 序列 的 小 毛 
ШЕ ЕО 

同 质 性 检验 结果 显示 已 >0. 01 ,说 明 125, 165 

































































统 发 育 树 中 , 弯 叶 毛 球 虫 属 和 食 螨 球 虫 属 则 是 舱 入 
到 小 毛 球 虫 属 里 面 (图 3) ,共同 形成 一 个 分 文 。 基 
于 联合 基因 构建 的 最 大 似 然 树 和 贝 叶 斯 推断 树 , 拓 
扑 结构 相似 (图 4) ,都 支持 了 小 毛 球 虫 属 为 单 系 起 
源 , 且 支持 率 较 高 (BS 296; PP =1) „(Нл EG 
内 各 亚 属 均 没 能 恢复 为 单 系 , 其 下 的 小 毛 际 虫 亚 属 、 
毛 球 虫 亚 属 . 小 球 虫 亚 属 和 拟 小 际 虫 亚 属 均 为 并 系 
起 源 。 这 与 基于 单 基因 12S, 165 和 28S rRNA 构建 
的 系统 发 育 树 相同 (图 1 ~4)。 系 统 发 育 树 主要 形 
成 6 KDEA), 16773 А 中 仪 包 含 小 又 虫 亚 属 ， 
支持 率 较 高 (BS =100; PP =1) , 另 有 部 分 小 球 虫 
亚 属 的 种 类 分 布 在 分 文 C, EMF, ÆI D 中 
仅 包 含 小 毛 际 虫 亚 属 , 但 仍 有 部 分 小 毛 际 虫 亚 属 的 
种 类 分 布 在 分 支 C, EA FP, ES HIE JE N 290 
分 布 在 分 支 B, E 和 下 中 ,也 没 能 形成 单 系 。 在 仅 涉 
及 的 拟 小 又 虫 亚 属 3 个 样本 中 , 立 拟 小 球 虫 S. 
(Pa. ) secula ^H Bp ШАХ SC o 5. (Ра.) 
alishanensis 聚 为 一 支 , 且 支 持 率 较 高 (BS =98; PP = 
1) ЖА ЛЕЗ B HP, -IUR JE AEIR 
亚 属 的 部 分 种 类 形成 一 个 分 支 , 另 一 个 样本 南 岭 拟 
Л 8 S. (Ра. ) nanlingicus Wii AIAX Е, 5 
KENELLE 5. (5. ) kabakovi R 73 — xc (BS = 
100; PP=1) 。 此 外 ,结果 还 表明 在 属 级 阶 元 的 系统 
发 育 关 系 上 , ESI UB Um 5 5 ESI Jy H OV 4H. 
妹 群 。 






















































































3 讨论 


线粒体 12S 和 16S rRNA 及 核糖 体 28S rRNA Ж 
因 是 昆虫 系统 学 研究 中 一 类 有 效 的 分 子 标记 , 因 其 
不 同 区 域 有 不 同 的 进化 速率 ,可 探讨 不 同 阶 元 的 系 
统 发 育 关 系 ( 朱 保 健 等 , 2010; 王 德 明 , 2010; 明 庆 

















和 28S rRNA 基因 之 间 系 统 学 信息 差异 不 显著 ,可 
进行 联合 序列 分 析 。 首 先 基于 单个 基因 分 别 构建 系 
统 发 育 树 ,利用 MrModeltest 3.7 АФ Be e 36 pi 
型 ,12S, 16S 和 28S rRNA 3 个 基因 均 选 用 GTR + 
I+G 模 型 用 于 最 大 似 然 树 和 贝 叶 斯 推断 树 的 构建 。 
然后 通过 PartionFinder 1. 1. 1 搜寻 联合 基因 最 佳 的 
分 区 策略 和 碱 基 替 换 模 型 ,结果 表明 每 个 基因 独立 
分 区 , GTR € I € G 为 每 个 分 区 最 佳 的 碱 基 蔡 换 
模型 。 

在 基于 12S rRNA 和 16S rRNA 基因 序列 构建 
的 系统 发 育 树 中 ES SCENES BS PR RE TR 
(图 1, 2) ,而 在 基于 28S rRNA 基因 序列 构建 的 系 


















































4845, 2013) 。 本 研究 依据 Magro 等 (2010 ) 37360 
IH 12S, 165 41285 rRNA 基因 的 引物 信息 ,成 功 扩 
В Л E SIL 9 Jg 125 rRNA 基因 约 356 bp、16S 
rRNA 基因 约 351 bp #1285 rRNA 基因 约 315 bp ,与 
Magro & (2010) 的 扩 增 结果 相 一 致 。 扩 增 的 线粒体 
125 rRNA #1 165 rRNA 基因 的 A+T 含 量 明 显 高 于 
G+C 含 量 , 表 现 出 明显 的 A +T 碱 基 偏 哮 性 ,而 且 
A 和 了 的 碱 基 含 量 相 当 , 符 合 Simon 报道 的 昆虫 碱 
ЖФА +T 比例 高 的 特点 (Simon et al., 1994) , 
Aruggoda 等 (2010) 通 过 扩 增 47 ЖЕГІҢ |Б) 165 rRNA 
基因 3 个 区 域 的 片段 也 表现 出 A +T 含 量 明显 高 于 
G+C 含量 的 特点 。 本 研究 中 核糖 体 基 因 285 rRNA 
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图 1 ЖТ 125 rRNA JEDE JJ EEM ERK СЕ BN 
Fig. 1 Phylogenetic tree of Scymnus based оп 125 rRNA gene sequence 


分 别 使 用 MrBayes 3.2.6 fll RAxML 8.2. 8 构建 的 系统 发 育 树 ,GTR +1 + G 作为 碱 基 最 佳 的 替换 模型 ;节点 旁边 前 面 为 自 举 值 , 后 面 为 后 验 概 
率 ; 图 2 和 3 F], The phylogenetic tree was reconstructed by MrBayes 3.2. 6 апа RAxML 8.2.8, and GTR +1 + С was used as the optimal substitution 
model. Bootstrap values and posterior probabilities for each branch are shown close to the nodes. The same for Figs. 2 and 3. 




















2 期 黄 卫 东 等 : 基于 12S, 165 和 28S rRNA 基因 序列 的 中 国 小 毛 融 虫 属 系统 发 育 分 析 








249 





92/1 


87/0.97 


73/0.84 


55/0.46 


23/0.42 


75/0.56 


80/0.57 


NE Ji Subgenus 
BEILER S. (5) 
园 小 球 虫 亚 属 S. (P.) 
Bes RS. (М) 

国 拟 小 又 虫 亚 属 S. (Ра) 
国 直 小 毛 际 虫 亚 属 9. (0.) 









































Fig. 2 
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Phylogenetic tree of Scymnus based on 16 rRNA gene sequence 
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基因 的 G+C 含量 略 高 于 A+T 含 量 , 这 与 已 有 的 研 
究 结果 相 一 致 ( 江 世 宏 等 , 2009; 张 健 等 , 2012; ІҢ 
REESE, 2013). 

在 基于 125 rRNA 和 16S rRNA 基因 序列 构建 
的 系统 发 育 树 中 , 均 表 明 小 毛 球 虫 属 为 单 系 起 源 


(图 1 和 2)。 然 而 ,在 基于 28S rRNA 基因 序列 构建 
的 系统 发 育 树 中 小 毛 标 虫 属 没 能 恢复 为 单 系 , 弯 叶 
EILE EAA HAUR JERA EDE BUE JE E ат, H 
同形 成 一 个 分 支 ( 图 3)。 可 能 是 由 于 28S rRNA Ж 
因 在 进化 过 程 中 比较 保守 ( 江 世 安 等 , 2009; 张 健 
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图 3 基于 28S rRNA ЫАЖ ERG ХАБ ІҢ 89 RRR BA 
Fig. 3 Phylogenetic tree of Seymnus based on 28S rRNA gene sequence 


等 , 2012) ,单独 作为 分 子 标记 用 于 际 虫 科 系统 发 育 系 。 然 而 ,28S rRNA 基因 与 其 他 分 子 标 记 联 合 使 用 
研究 时 ,未 能 反映 出 际 虫 科 属 级 阶 元 的 系统 发 育 关 ”可 以 较 好 地 解决 昧 虫 科 不 同 阶 元 的 系统 发 育 关 
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Fig. 4 Phylogenetic tree of Scymnus based on the concatenated data of 125, 165 and 28S rRNA gene sequences 
分 别 使 用 MrBayes 3.2.6 ЖІ RAxML 8. 2. 8 构建 的 系统 发 育 树 ,每 个 基因 独立 分 区 ,GTR +1 +G ТЕУ ІК БОНН ФАН; ;节点 旁边 前 面 为 


















































自 举 值 ,后 面 为 后 验 概率 。The phylogenetic tree was reconstructed by MrBayes 3. 2. 6 and RAxML 8.2.8. Each gene was zoned separately апа GTR + 


I +G was used as the optimal substitution model for the three gene partitions. Bootstrap values and posterior probabilities for each branch are shown close 


to the nodes. 
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系 (Giorgi et al., 2009; Magro et al., 2010; Seago et 
al., 2011; Szawaryn et al., 2015) 。 

联合 基因 构建 的 最 大 似 然 树 和 贝 叶 斯 推断 树 都 
支持 小 毛 驹 虫 属 为 单 系 起 源 (图 4) ,这 与 Giorgi 等 
(2009) „Margo 等 (2010)、Aruggoda 5 (2010) „Seago 
^ (2011) 和 Robertson 45 (2015) 的 研究 结果 相 一 
$t, [HE TERI АЙЯ ҰНЫ, (УР RIP EE R^ NE SUR 
属 样本 ,代表 性 不 强 ,而 本 研究 首次 基于 广泛 的 取样 
证 实 了 小 毛 际 虫 属 为 单 系 起 源 。 另 外 ,在 小 毛 标 虫 
属 内 ,5 个 亚 属 之 间 的 系统 发 育 关系 较为 混乱 。 直 
小 毛 飘 虫 亚 属 仅 选择 了 1 种 ,无 法 讨论 该 亚 属 的 单 
系 性 ,而 其 他 4 个 亚 属 均 未 恢复 为 单 系 。 此 外 ,在 关 
于 球 虫 科 的 系统 发 育 研 究 中 ,多 数 研究 者 对 于 小 毛 
球 虫 属 并 没有 采用 亚 属 的 分 类 系统 ,只 有 Aruggoda 
等 (2010 ) 的 研究 中 采用 亚 属 的 分 类 系统 , 仅 涉 及 2 
个 样本 :肥厚 小 毛 早 虫 S (S. ) pinguis AUS Д 
H S. (Р. ) oestocraerus ,无 法 阐明 各 亚 属 之 间 的 系 
统 发 育 关系 。Slipifiski 等 (2012 ) 系统 研究 了 澳 大 利 
亚 的 小 毛 球 虫 属 , 认 为 目前 该 属 的 亚 属 分 类 特征 与 
澳大利亚 的 小 毛 标 虫 种 类 不 相符 ,因此 舍弃 各 亚 属 
而 采用 了 广义 上 的 小 毛 标 虫 属 。 

在 小 毛 飘 虫 属 传统 的 形态 学 研究 中 ,主要 依据 
触角 节 数 .后 基线 是 否 完整 以 及 弯 管 圳 的 形态 结构 
为 主要 特征 来 区 分 这 5 个 亚 属 ( 陈 晓 胜 , 2013 ) 。 本 
研究 的 结果 显示 上 述 5 个 亚 属 均 未 能 恢复 为 单 系 ， 
这 表明 小 毛 际 虫 属 下 面 的 亚 属 分 类 系统 需要 重新 修 
订 。 同 时 ,也 表明 仅 通过 上 述 形态 指标 用 于 区 分 亚 
属 仍 有 较 大 的 局 限 性 ,因此 ,在 比较 形态 学 的 研究 
中 ,应 该 挖掘 更 多 的 形态 特征 ,重新 构建 各 亚 属 的 形 
态 指标 体系 。 

关于 小 毛 球 虫 属 S EIU E RU Er А CB, Ji 
三 者 之 间 7 NES 5 2$ m E91 8 Ја [8] 88 FE 
ЕЛЕР , m AWAUE JE DU] Je T UR S TENE 
AEDSTA EE RU SS В, RUE] TEAU S 5 
St E AUR ЖАС ЖОЛ ( Magro et al., 
2010) , MA MAUR JE 55 25 n -EL LB + АД 
属 互 为 姐妹 群 (图 4) ,这 与 基于 形态 学 数据 的 研究 
结果 相 一 致 ( Yu, 1994); 

本 研究 首次 基于 线粒体 基因 128 和 165 rRNA 
及 核糖 体 28S rRNA 基因 共 3 个 分 子 标记 ,使 用 ML 
法 和 BI 法 重建 了 我 国 小 毛 标 虫 属 亚 属 的 系统 发 育 
关系 ,但 其 下 的 亚 属 多 为 并 系 起 源 , 且 亚 属 之 间 杀 缘 
关系 较为 混乱 。 由 于 物种 有 限 ,并 不 能 精确 地 解析 
小 毛 飘 虫 属 内 各 亚 属 之 间 的 系统 发 育 关 系 。 因 此 关 










































































































































































于 小 毛 疆 虫 属 内 各 亚 属 的 分 类 地 位 是 否 成 立 ,还 需 
要 做 进一步 研究 。 但 本 研究 测定 的 44 种 小 毛 昧 虫 
HS 125, 165 和 28S rRNA 基因 序列 为 后 续 小 毛 又 
虫 属 的 进一步 系统 发 育 研 究 芮 定 了 一 定 基 础 。 
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